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Abstract 

We present the one-loop helicity amplitudes with five external gluons. The computation em- 
ploys string-based methods, new techniques for performing tensor integrals, and improvements in 
the spinor helicity method. 
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Calculations beyond the leading order in quantum chromo dynamics are important to refining 
our understanding of known physics in present-day and future collider experiments, such as the 
Tevatron or the SSC and LHC. In jet physics, next-to-leading order calculations are important 
to curing several deficiencies of their leading-order counterparts: the strong spurious dependence 
on the renormalization scale; the lack of sensitivity to the jet resolution parameters, namely the 
minimum transverse energy and the jet cone size; and the absence of a warning about dangerous 
infrared logarithms. The one-loop corrections to matrix elements for 2^2 processes in QCD, a key 
ingredient of the next-to-leading order calculations of inclusive-jet and two-jet cross-sections and 
distributions, were computed by Ellis and Sexton [1]. To go beyond these cross-sections, whether to 
higher orders for two-jet cross-sections, or to next-to-leading order for three-jet cross-sections and 
distributions, requires the computation of the one-loop corrections to the 2^3 matrix elements. 
At hadron colliders, the QCD coupling a s , and the manner of its running, can be extracted from 
purely hadronic processes by comparing three-jet production to two-jet production, at various 
center-of-mass energies. The presence of infrared logarithms in both of these quantities means that 
this cannot be done sensibly unless both quantities are known to next-to-leading order. 

We present here the one-loop matrix elements with five external gluons, which are the hardest 
part of a 2 — > 3 calculation if a traditional diagrammatic method is used. We have performed the 
calculation using the string-based methods developed in ref. [2] as more efficient tools for one-loop 
calculations with external gluons. The rules presented there were derived by taking the infinite- 
tension limit of an appropriately-constructed heterotic string amplitude. The structure of the rules 
can also be understood in conventional field theory [3], and the application of such methods to 
a calculation such as the present one does not require any knowledge of string theory. (It turns 
out that it is possible to construct a set of rules yielding more compact integral representation of 
gluon amplitudes at intermediate stages than would emerge from a straightforward application of 
the rules in ref. [2]. This alternate set of string-based rules will be discussed elsewhere.) 

In the string-based method, one first decomposes the n-gluon amplitude, depending on the 
external momenta, helicities, and color indices fcj, Aj, and m, into sums over certain permutations 
of color factors, times partial amplitudes, in analogy to the helicity [4,5] and color [6] decomposition 
of tree amplitudes. At one-loop order in an SU (N) theory, one must also sum over the different 
spins J of the internal particles; this takes the following form when all internal particles transform 
as color adjoints, 

_ |n/2J+l _ 

A n {{k\ i ,a i }) = Y j nj E Gr n;c (a) A^ c (a) (1) 

J c=l a€S n /S n . c 

where Gr n; i(l) = N Tv (T ai ■ ■ -T a "), Gr n;c (l) = Tr (T ai ■ ■ ■ T^- 1 ) Tr (T a - ■ ■ • T a " ) , S n is the set 
of all permutations of n objects, and S n;c is the subset leaving the trace structure Gr n;c invariant. 
The T a are the set of hermitian traceless N x N matrices, normalized so that Tr (T a T b ) = S ab . For 
internal particles in the fundamental (N + N) representation, only the single-trace color structure 
(c = 1) is present, and it is smaller by a factor of N. We take in each case a spin- J particle with 
two states: gauge bosons, Weyl fermions, and complex scalars. 

The objects one calculates are the partial amplitudes An-} c , which depend only on the external 
momenta and helicities. For the five-point function, there is only one independent partial amplitude 
for each configuration of external helicities; A 5;2 and A 5 . 3 are related to the adjoint contributions 
to Az, ; i via decoupling equations [7]. 

The string-based method meshes well with the spinor helicity representation for the polariza- 
tion vectors [4,5], which provides an efficient method for extracting the essential gauge-invariant 
information in an on-shell amplitude. This method yields expressions written in terms of spinor 
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products (ij) and [ij], which are square roots of Lorentz products Sij = (ki + kj) 2 (up to a phase). 
Unfortunately, the relations — momentum conservation and the Schouten identity — between dif- 
ferent forms of a given expression are nonlinear, which makes it hard to give a canonical form for 
such expressions, or equivalently makes it hard to simplify complicated expressions. However, one 
can evaluate "phase-invariant" combinations of spinor products in terms of Sij and contractions of 
the Levi-Civita tensor e(i,j, m, n) = 4ze M!/p(T /cf kjk^k^ = [ij] (j m) [mn] (ni) — (ij) \j m] (mn) [ni]. 
It suffices to calculate the ratios 

(ij) \j k] sgSjk + SkiSjj - s ik Sji - e(i,j,k,l) 

(it) [Ik] ~ 2s u s kl ' 1 ' 

using e.g. methods in reference [8]. In this way spinor products can be eliminated from any 
expression, apart from an overall prefactor. 

For massless five-point kinematics, such an expression can then be written as a rational function 
in the five kinematic variables {Pi, P3, (3%, P$} (or any cyclic permutation of this set), where 

5 1/2 

Pi = [i, i + + » + 2)[i + 2, i + 3](i + 3, i) x (- JJ*ij+i) • ( 3 ) 

The only independent Levi-Civita contraction is given by e(l, 2, 3,4)/(— n*=i s j,j+i) 1 ^ 2 = (0502 + 
P1P3 + P2P4 + P3P5 + P%Pi)IPz = Pi — Pi for any i, and the independent Lorentz products by 
Sj 5 j + i = —l/((Pi + P* +1 )(Pi + 2 + P* + z))- Simplification of rational functions in Pi is straightforward. 

The Pi variables are related to the variables ji and A 5 used in ref. [9] to perform pentagon 

integrals, via P^ = — (7^+2 ± V^5)/2. Indeed, the derivative approach to evaluating tensor 
integrals [9], when applied to the pentagon integrands encountered in the five-gluon calculation, 
and expressed in terms of the appropriate set of Pi variables, allows one to significantly reduce the 
degree and size of the Feynman parameter polynomials in the integrand. 

At tree- level, certain helicity amplitudes vanish identically [10]. The corresponding one-loop 
amplitudes are then free of infrared divergences. The remaining amplitudes are infrared-divergent; 
for practical purposes these divergences must be regulated using dimensional regularization. The 
computation of these helicity amplitudes thus requires the knowledge of five-point loop integrals in 
D = 4-2e [9, 11]. 

For the finite helicity amplitudes, supersymmetric identities [12] imply that the contributions 

of particles of different spin circulating around the loop are related, A$ c = -A [ l{? ] = A [ °) c . (This 
holds true for the partial amplitudes whether or not the theory as a whole is supersymmetric.) 
Indeed, in the string-based method, these identities hold for the integrands of each diagram. The 
amplitudes are 

/ + + + .+ +x _ i S12S23 + S23S34 + S34S45 + S45S51 + S51S12 + e(L 2, 3, 4) 
5:1 [ ' ' ' ' ' ~ 96vr 2 (12) (2 3) (3 4) (4 5) (51) 



4 1 1 (l-,2+,3 + ,4+,5+) 



(S23 + «34 + s 45 ) I 2 5 ] 2 - [ 2 4 ] (4 3) [3 5] [2 5] 



48tt 2 [12] (2 3) (34) (4 5) [51] 

_ Mj4 ( (12) 2 {13) 2 M- + (13) 2 (14) 2 14 + (14> 2 {15) 2 M.) 
(12) (15) \ (23) X 1 V 1 (34) x ' x 1 (45) J 

(4) 

In order to present the results for the remaining, infrared-divergent amplitudes in a compact 
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form, it is helpful to define the following functions, 



Lo(r) 



Lsi(ri,r 2 ) 



ln(r) 
1^7' 



Li(r) 



ln(r) + 1 — r 
(1 -r) 2 



L 2 (r) 



ln(r) - (r - l/r)/2 



Li 2 (l - ri) + Li 2 (l - r 2 ) + lnri lnr 2 - 



7T 



(1 - n - r 2 ) 2 

+ (l-ri-r 2 )(Lo(ri)+Lo(r 2 )) 



(5) 



where Li 2 is the dilogarithm; a prefactor, 



cr 



(47r) e r(l + e)r 2 (l-e) 
16vr 2 T(l - 2e) 



(6) 



a universal function, 
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' ^ ' + £ In ( -=^) In ( ^±^i) + V - ^ , (7) 



the following functions for the (1 ,2 , 3 + , 4 + , 5 + ) helicity configuration, 



vf 

F f 



5 1 
Ye ~ 2 



In 



-S23 



+ In 



-S51 



+ - 
3 9 



-S51 



1 (1 2) 2 ((2 3) [3 4] (4 1} + (2 4} [4 5] (5 1)) L o 

2 (2 3) (34) (4 5) (5 1) 

1 [34] (41) (2 4) [4 5] ((2 3) [34] (41) + (2 4) [4 5] (51)) L2 (fff j | 



S51 



(3 4) (4 5) 



51 



(3 5) [3 5r 



+ 



1 (12) [3 5]' 



+ 



1 (12) [34] (41) (24) [4 5] 



3 [12] [2 3] (3 4) (4 5) [51] 3 [2 3] (3 4) (4 5) [5 1] 6 s 23 (3 4) (4 5) s 5 i 



(8) 
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and the corresponding ones for the (1 , 2 + , 3 , 4 + , 5 + ) helicity configuration, 



5_ 1 

Ye ~ 2 



In 



-S34 







+ ln| 


' ^ 






-«23 


-S34 ] 


— 851 ' 


-S51 / 


s 2 
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1 , 2 

-2. v = --f' + - 



(13) 2 <41}[24] 2 La > (3ff. 5Sf) <13} 2 <53}[25] 2 Ls i (^ 3S) 



F s = — 



(4 5} (51) si, (3 4) (4 5} S 2 4 

1 (13) 3 ((15) [5 2] (2 3) - (34) [4 2] (21)) Lo (gg) 

2 (12) (2 3) (3 4) (4 5) (51) a 51 

(12)(23)(34)(41) 2 [24] 2 2 ^1 (gg, gg)+Li (gg)+Li (gg) 
(4 5)(51)(24) 2 "~~ " " ~ s§i 

o\2r oc -!2 2 LSi 



(3 2) (21) (15) (5 3) 12 5 

"l" \2 



+ 



+ 



(54)(43)(25) z ^ 4 (g) 

2 (2 3) 2 (41) 3 [2 4] 3 L2(5fg) _ 2 (2 1) 2 (5 3) 3 [2 5] 3 L 2 (gg) 
3 (4 5) (51) (2 4) ^ 3 (5 4) (4 3) (2 5) sf 4 

L2 (5fff) /I (13) [24] [2 5] ((15) [5 2] (23) - (34) [42] (21)) 



4i V3 (4 5) 



2 (12) 2 (3 4) 2 (41) [2 4] 3 2 (3 2) 2 (1 5) 2 (5 3) [2 5] ; 



3 (45) (51) (24) 3 (54) (43) (25) 

| 1 (13) 3 ((15) [52] (23) -(34) [42] (21)) Mgg) < 1 [ 24 ] 2 [25] 2 



6 (12) (2 3) (34) (4 5) (51) s 51 3 [1 2] [2 3] [3 4] (4 5) [5 1] 



(12)(41)"[2 4] J 1 (3 2)(5 3H2 5r | 1 (1 3)' [2 4] [2 5] 



3 (4 5) (51) (2 4) [2 3] [3 4] s 51 3 (5 4) (4 3) (2 5) [2 1] [1 5] s 3i 6 s 34 (4 5) S 5i ' 

For positive values of Sij, the logarithms and dilogarithms develop imaginary parts according to the 
prescription — > Sij+ie. We also remind the reader of the tree amplitudes, A^ ee (1~, 2~, 3 + , 4 + , 5 + ) 
i(12) 4 /((l 2) (2 3) (34) (45) (51)) and A^ e (1", 2+, 3", 4+, 5+) = i{\ 3) 4 /((l 2) (2 3) (34) (45) (51)). 
In terms of these functions, the MS renormalized amplitudes are 

A™ = Cr (V s Af ee + iF s ) , 
A% 2] = -c r {(V f + V s )AT e + i(F f + F s )) , (10) 
= c r ((V 9 + 4T/ / + y s )yl| ree + i(4F f + F s )) . 

The rest of the helicity amplitudes are related by cyclic permutations or complex conjugation to 
those given above. It is interesting to note that in supersymmetric theories, the V s and F s terms 
cancel out of the final amplitude, and that in TV = 4 supersymmetric theories only the V 9 term 
survives. The separation implied above into g, /, and s pieces arises naturally on a diagram-by- 
diagram basis within the string-based approach. In this approach the V 9 term represents the only 
calculational difference between the contributions with gluons circulating around the loop, and 
those with fermions; this term has a particularly simple expression at every intermediate stage of 
the calculation. The parameter Sr controls the variant of dimensional regularization scheme [2]: 
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for 8r = 0, one obtains the four-dimensional helicity scheme, while for Sr = 1 one obtains the 
't Hooft-Veltman scheme. 

There are several checks we have applied. We have checked gauge invariance, both by com- 
puting amplitudes with longitudinal gluons, verifying that one obtains zero, and by calculating 
a helicity amplitude with an alternate choice of spinor-helicity reference momenta, and verifying 
that the result is unchanged. In addition, the forms given above display manifestly the reflection 
symmetries expected of the amplitudes, symmetries that are not present in the contributions of 
the individual diagrams. The amplitudes also have consistent limits as one of the gluon momenta 
becomes soft, and as two adjacent momenta become collinear. 

At next-to-leading order, only the infrared-divergent helicity amplitudes (5-10) enter into the 
construction of a program for physical quantities. In order to construct such a program for three-jet 
quantities, one must form the interference of the tree amplitude with the loop amplitude; this has 
the form [7] 



colors 



Re 4 ree >M5;i(^) 
<res 5 /z 5 



E I^AInlo = (iV 2 - 1) 

+ ^ Re E [A t r*(r.p)A 5;1 (p)-A t r*(p)A 5;1 (r-p) 



pes 5 /z 5 



(11) 



+^ Re E E Al-*(h-p)A 5;3 (p) 



heH 5 



where r is the permutation (2413 5); ^(3) is the ten-element set of distinct partitions of five ele- 
ments into lists of length two and three; and H 5 = {(12345), (3 4125), (3 1245), (21345), (3 2145), 
(34215)}. For QCD with rif flavors of massless quarks, one substitutes A^.\ — > A^\ + ^^^i 2 ' 

and A^- 3 — > A^) 3 into equation (11). One must then combine this virtual correction with the sin- 
gular terms in the 2^4 matrix elements arising from the integration over soft and collinear phase 
space. The Giele-Glover formalism [13] makes use of the color ordering in construction of universal 
functions representing the results of the soft and collinear integrations, and is the most convenient 
one for doing so. We have used it to check that the poles in e do cancel as expected. 
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